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Abstract 

Subcritical transition of an inhomogeneous plasma where turbulences with different characteristic 
space-time scales coexist is analyzed with methods of statistical physics of turbulences. We derived 
the development equations of the probability density function (PDF) of the spectrum amplitudes 
of the fluctuating electro-static potential. By numerically solving the equations, the steady state 
PDFs were obtained. Although the subcritical transition is observed when the turbulent fluctua- 
tions are ignored, the PDF shows that the transition is smeared out by the turbulent fluctuations. 
It means that the approximation ignoring the turbulent fluctuations employed by traditional tran- 
sition theories could overestimate the range where hysteresis is observed and statistical analyses 
are inevitably needed. 
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I. INTRODUCTION 



Transition phenomena with sudden changes of states are observed in turbulent plasmas. 
Since these transition phenomena like the L-H transition play crucial roles for magnetic 
confinement of fusion plasmas, the transition associated with formation of transport barriers 
is one of the main subjects of high-temperature plasma physics. 

Traditional theories for transition of high-temperature plasmas are formulated in terms 
of averaged physical quantities Fluctuations around the averages are ignored and 

transition phenomena are described deterministically. 

However, high-temperature plasmas are strongly non-linear systems with huge number 
of degrees of freedom and hence their behavior should be chaotic and unpredictable in a 
deterministic way. In fact, broad distribution of critical values of parameters where transi- 
tions occur and intermittent transport called "avalanche phenomena" are observed in recent 
experiments 3, U, 0]. Occurrence of these behaviors, which cannot be described only by 
averaged quantities, are considered due to strong turbulent fluctuations. It is inevitable to 
describe turbulent plasmas in terms of probability or ensembles like numerical forecast of 
weather, since the magnitudes of fluctuations around averages are of the same order as that 
of the averages in turbulent states. 

Therefore, transition takes place as a statistical process in the presence of stochastic noise 
sources induced by turbulence interactions. As a generic feature, transition is expected to 
occur with a finite probability when a control parameter approaches the critical value. 

Statistical theories for plasma turbulence have been developed and the framework to 
calculate the probability density function (PDF), the transition probability etc. has been 



made 
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121 ] . In the statistical theories, the time-development of the system 



is described by a set of differential equations with random forces, called the "Langevin 
equations". All the information on the statistical properties of the system is obtained by 
solving the Langevin equation. 

The framework has been applied to cases where only one turbulent mode is excited and 
the turbulence is characterized by one space-time scale jy, |3, 

IU 114113- 

However, it is well known that there are many kinds of turbulent fluctuations in high- 
temperature plasmas and that different characteristic length scales coexist. The importance 
of interactions between modes with different scale lengths has recently been recognized. For 
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instance, the dynamics of the meso-scale structure of the radial electric field [16, [17|, ll8| is 



known to cause variation in the dynamics of microscopic fluctuations like in the electric filed 



domain interface 



19, 201, zonal flow 



2l| and streamer 



221 ] . Coexistence of multiple sea 
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turbulence has also been investigated by use of the direct numerical simulations 
Statistical theory on zonal flow dynamics j3| and that of the L-H transition theory has been 
developed 3|. 

In the present paper, we apply the statistical theoretic algorithm to a model of high- 
temperature plasma where two characteristic scales coexist. One is the current diffusive 
interchange mode (CDIM) micro turbulence Q|, whose characteristic length scale is of the 
order of the collisionless skin depth 5 = c/uj p . The other is the ion-temperature gradient 
(ITG) mode turbulence, whose characteristic wave length is of the order of the ion gyroradius 
pi, as an example of the drift wave fluctuations considered to dictate a considerable part 
of the turbulent transport [26] . Hereafter, we call these two modes "the high wave number 
mode" and "the low wave number mode" respectively for its simplicity. We assume that the 
condition pi 3> 5 holds. Both turbulences are considered to cause the anomalous transport 
and hence the coexistence of the high wave number mode and the low wave number mode 
turbulences and their interplay should be taken into account. 

It is known that the subcritical transition occurs in this system, when the pressure- 
gradient and the radial electric field are changed j27(. However, the turbulent fluctuations 
are ignored in the analysis. In the present paper, with the statistical theory, we analyze 
the effect of the turbulent fluctuations on stochastic properties of the transition and show 
that the fluctuation changes the phase structure of the system completely. More precisely, 
we show that the transition is smeared out by noises, i.e., the physical quantities changes 
gradually without clear transition. 

The present paper is organized as follows; the statistical theory and the model are formu- 
lated in Sec. |H] The results of deterministic analyses including occurrence of the subcritical 
bifurcation are summarized in Sec. IIHI The Section HVj i.e., the main part of this paper, 
presents the statistical properties of the system. Summary and discussions are given in Sec. 
lYl 
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II. THEORETICAL FRAMEWORK AND THE MODEL 



In this section, we briefly review the theoretical framework and the model of turbulent 
plasmas where two different characteristic scales coexist 27 1. 

The starting point of the theory is the reduced MHD for the three fields: the electro-static 
potential, the current and the pressure. The Langevin equation that the statistical theory is 
based upon is derived as one of model equations which reproduce the two-time correlation 
functions and the response functions obtained by the renormalization perturbation theory 
(the direct-interaction approximation) for the reduced MHD j^. 

The Langevin equation describes time-development of two variables characterizing the 
system. One is the spectrum amplitude of the electro-static potential for the characteristic 
wave number kh of the high wave number mode, Xh = k\{(f)* kh (j)k h ) / 'D^, and the other is that 
for the low wave number mode, Xi = kf (0£ 0*;,) / Df , whose characteristic wave number is k\. 
Here, and Di denote the renormalized transport coefficients for the high wave number 
mode and the low wave number mode respectively when there is no interactions between 
two modes. The characteristic time constants for the two modes are defined as 

ut = k?D h u h = k 2 h D h . (1) 

See [27] for the details of the notation. 

The Langevin equation gives the time-development of these two variables as 



d(ujit) 2 

dx h 



dxi 1 v /x^ jx h \ ~ n \ __ If y/x^ 



d{u h t) 



+ 



i% + - 2 )^ + 2\r + \l^ + 4Xl ) XlWl = °' (2) 

x h + x 3 h /2 w h + y/ewt = 0. (3) 



1 + /(l + pr 2 Xl 



\ y/x~h/r + \Jx h /r 2 + Ax{ 
Nonlinear drags and magnitudes of the noises have been evaluated in 

In these equations, the nonlinear terms of the reduced MHD are divided into two parts: 
One part is coherent with the test field and is renormalized into the deterministic terms, 
i.e., the second terms. The other is incoherent and is modeled by random forces Wh,wi. 
Another random force Wt denotes the thermal noise and the magnitude of the noise, e, is a 
small quantity compared to the magnitudes of other noises. Technically, introduction of e is 
needed in order to exclude singularity at Xh = of the PDF. 

Since Wh,wi and w t are forces which fluctuates randomly in time, the Langevin equation 
describes the stochastic time-development of the turbulent fluctuation of the system. For 



simplicity, we assume that the random forces are Gaussian and white; 

(w(t)w(t')) = 5(t-t>). (4) 

This system is sustained by the space inhomogeneities, the curvature of the magnetic field 
Q', the pressure gradient dPo/dx and the gradient of the radial electric field dE r /dr. Here, 
the shear of the magnetic field is given as the slab configuration: B = (0, B sx, B ) where 
Bq{x) = const x (1+Q'x+- ■ •). The pressure is assumed to change in x— direction. These driv- 
ing forces are characterized by the parameter r and p as r = Di/Dh oc {Vt'dP^/dx) -1 ^ 2 ,p = 
Dh/Ieff 1 oc (dE r /dr)~ 6 Q'dPo/dx. Here, I^ft denotes the critical strength of the nonlinear 
interactions between the low and high wave number modes. The pressure gradient controls 
the growth rate of the low wave number mode and Q'dPo/dx excites both the high and 
low wave number mode turbulences. The gradient of the radial electric field suppresses 



turbulences 
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We assume that the relation 3> 5 holds and hence the characteristic length-scales for 
the two modes are widely separated as 

k h > k. (5) 

The mutual interactions between the low and high wave number modes are asymmetric, 
since the spatial structure of the low wave number mode is a large-scale inhomogeneity for 
the high wave number mode. The assertion Eq. (JSJ also means that the time-scales are 
widely separated since the time scales are given by Eq. (JTJ); 

u h > uji. (6) 

By analyzing the Langevin equation, Eqs. (0 El), a number of statistical properties of 
turbulent plasmas can be derived. For example, the analytical formulae of the rate of change 
of states of plasmas, the transition rates, were derived. Furthermore, since the renormalized 
transport coefficients and the random forces have the same origin, i.e., nonlinear interactions 
in MHD turbulence, relations between the fluctuation levels of turbulence and the transport 
coefficients like the viscosity and the diffusivity were derived. 



III. BIFURCATION WITHOUT RANDOM FORCES 



With the theoretical framework briefly described in the previous section, we analyze the 
model of the inhomogeneous plasma with two characteristic scales. 
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FIG. 1: The r-dependence of the steady state solutions when no random force. The parameter p 
is fixed as p = 30. The solid line represents the low wave number mode x\ and the dotted line 
represents the high wave number mode Xh- It is seen that the subcritical bifurcation occurs. 

At first, we show the steady state solutions when random forces are ignored, in order to 
compare to the results with the random forces obtained later. The steady state solutions 
are obtained by solving the set of nonlinear equations when the random forces are turned 



Figure H shows the r-dependence when p = 30. When r < 0.49 . . ., the low wave number 
mode turbulence is suppressed. As r is increased, the system experiences the subcritical 
transition to the state where the low wave number mode turbulence is excited. When 
0.49 . . . < r < 1, there are two stable solutions and it means that the system is bi-stable. 
From the deterministic point of view, the transition is expected to occur only at the ridge 
point and the bifurcation point. The qualitative behavior does not depend on the value of 
p as far as p > 1.9 . . .. 

The phase diagram is shown in Fig. 121 The subcritical bifurcation is observed when the 
value of the parameter p is larger than 1.9 . . .. On the other hand, when p < 1.9 . . ., the 
bifurcation is supercritical. 
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/(l+pr 2 xi) x h = 0. 
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FIG. 2: The phase diagram when the random forces are turned off. The subcritical bifurcation 
is observed when p > 1.9... and the boundaries of the bi-stable phase are represented with the 
solid line. When p < 1.9 . . ., the bifurcation is supercritical and critical values are plotted with the 
dotted line. The region indicated as "h phase" is the state where the high wave number mode is 
excited and the region "1 phase" is the state where the low wave number mode is excited. 

IV. THE STOCHASTIC PROPERTIES 

In the rest of the present paper, we analyze the stochastic properties of the model, Eqs. 
(jSHSJ)? to investigate the effect of the turbulent fluctuations, i.e., the random forces. 

A. The adiabatic approximation and the Fokker-Planck equation 

At first, we approximate the Langevin equation Eqs. (J21 EJ) with making use of the time 
scale separation between the low and high wave number mode turbulences. 

The scale separation, Eq. means that the high wave number mode Xh quickly relaxes 
to the steady state determined by the value of the low wave number mode variable x\ which 
is fixed at the value xi = xi(t) at the time. 

We analyze the steady state of the high wave number mode Xh when xi is fixed. A state 
of a stochastic system is described by the probability that a state variable takes a certain 
value. The time-development of the probability density function (PDF) P(xh,t) of Xh is 
determined by Kramers-Moyal expansion applied to the Langevin equation of Xh, Eq. Q. 
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The Kramers- Moyal expansion is given as 
dP(x h ,t) d 



1 d 2 



C x {x h )P{x h ,t) + --^jC 2 (x h )P(x h ,t). 



(9) 



dt dxh 

The expansion can be truncated at the second order, since the random forces in the Langevin 
equation are assumed to be Gaussian. Here, the coefficient C n (xh) is given by 



C n (x h ) = hm i ([x h (t + At) - x h (t)] n ) 



(10) 



where (■) denotes the average over the all realizations of the random forces and the average 
is taken under the condition Xh{t) = x^. The resulting equation of motion of the probability 
density function is called the Fokker-Planck equation and is written as 



dP(x h ,t) _d_ 
d{u h t) dx h 



x h /r+ sjx h /r 2 + Axi 

i d 2 



/(l +pr 2 X[) 



xx h P(x h} t) + - ^-(4 + e)P( x h, t). 



2dx 2 



(11) 



The steady state solution of Eq. (JTTj) when x\ is fixed is shown in Fig. El It is seen that 
the peak of the PDF is relatively narrow and it means that the high wave number mode 
turbulence spends most of time at the peak value x* h . Hence, we can say that it is good 
approximation to replace Xh in the Langevin equation of the low wave number mode with 
the peak value x* h . 

The equation to determine the location of the peak x* h is given by the condition 

dP(x* h )/dx* h = as 



\ 



i + 



' X l If -i i 2 \ 3 

^_ / ; = /(I + P r x ~ o X h ~ 

r x\jr + Jx* h /r 2 + Ax i A 



0. 



(12) 



It is important to note that Eq. ()12j) is essentially different from Eq. (jHJ) in existence of the 
second term of Eq. (|12jh which comes from the random force x^Wh of Eq. Q. It implies 
that the random forces change the steady state of the high wave number mode. 

Consequently, the adiabatically approximated Langevin equation for the low wave number 
mode is given by Eq. (j2J) where Xh is replaced with x* h determined by Eq. (JTJJ). The reduced 
Langevin equation is written as 



dxi 1 
d(u t t) + 2 



'xl 



+ 



+ Axi 




+ Axi xiwi = 0. 



(13) 



16 



14 - 



Q. 




12 - 



10 - 



2 - 



4 - 







6 - 



i 








0.2 0.4 0.6 0.8 



1.2 



1.4 



FIG. 3: The steady state PDF P(xh) when xi is fixed at x\ = 1. The peak of the PDF is relatively 
narrow and the adiabatic approximation that Xh is replaced with the peak value of the PDF, x^(xi), 
is justified. The arrow indicates the location of the peak. 

The corresponding Fokker-Planck equation, which determines the time-development of the 
PDF of xi, P(xi,t), is given by 



B. The probability density functions and the effect of the random forces 

In this subsection, we investigate properties of the steady state PDF of the low wave 
number mode, P(xi), by numerically solving Eq. (jl4j) . Although we will show figures when 
the parameter p is fixed at p = 30, qualitative behavior is the same as that for other p values 
larger than 1.9 . . .. 

At first, we show in Fig. 01 the steady state PDF in the small r region. In this region, 
we have seen that the low wave number mode is suppressed when the turbulent fluctuations 
are ignored. Figure 0] shows that the probability that the low wave number mode is quiet is 
large and hence the analysis ignoring random forces is a good approximation in this region. 

However, as the value of the parameter r is increased, the characteristics change com- 




(14) 
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FIG. 4: The steady state PDF P{x{) in the small r region. It shows that the probability that 
the low wave number mode is suppressed is large and the result is compatible with that of the 
deterministic analyses given in Fig. ^ 

pletely. The steady state PDF when r > 0.49 ... is shown in Fig. El In the region 
0.49... < r < 1, we have seen that the subcritical transition occurs and the system is 
bi-stable if the random forces are ignored. Although there are two peaks in the PDF which 
are compatible with the previous result without random forces, the valley between the peaks 
is too shallow to identify the two states. In other words, even if one observes the time-series 
of xi(t), the value of xi(t) strolls around the two peaks without a sudden change. It means 
that the bi-stability of the system, i.e., the subcritical transition, is smeared out by the 
turbulent fluctuations. 

The steady state PDF for r > 1 is also shown in Fig. It is seen that the PDF has 
a single peak. The single peak is compatible with the result of the deterministic analysis, 
where only one state exists. However, the peak is wide and hence xi fluctuates widely around 
the location of the peak. The variance of the fluctuation is as large as the peak value and it 
means that the fluctuation cannot be ignored. Furthermore, the peak value 0.29 . . . obtained 
when r = 1.5, p = 30 is different from the root of Eqs. ([ZJIBJ), which is given as x\ ~ 0.98. 

In order to see these consequences explicitly, in Fig. we show the comparison of the 
contour plot of the PDF and the bifurcation diagram obtained by the analysis ignoring the 
random forces. It is seen that the location of the region which has large probability and 
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x, 

FIG. 5: The steady state PDFs of x\ when 0.49. . . < r. The solid line represents the PDF for 
r = 0.7 and the dotted line represents the PDF for r = 1.5. The wide peak which consists of 
two small unresolved peaks means that the subcritical transition cannot be observed due to the 
turbulent fluctuations. 

mainly denoted with the thin solid lines in Fig. |H1 is shifted gradually as the value of the 
parameter r is changed. Furthermore, the location of the peak of the PDF is far from the 
result when the random forces are ignored when r > 0.49 .... 

Next, in order to investigate the meaning of the average in such a turbulent system, we 
analyze the tail of the steady state PDF P(xi) shown in Fig. It is shown that the tail 
of the PDF is well-approximated by the power-law and hence the probability is distributed 
broadly over large values of the fluctuation. The exponent is about —0.34 for the case shown 
in Fig. [7| It means that the average, i.e., "the center of mass" of the PDF, is shifted in the 
large xi direction with the long tail and hence the average is not equal to the value that x\ 
takes with large probability. The values that xi takes with large probability are characterized 
by the peaks of the PDF and are called the most probable values. The difference between 
the average and the most probable value is shown in Fig. |H| 

Furthermore, hysteresis cannot be captured by observing the average, since the quantity 
is single-valued from its definition. On the other hand, the most probable values depict the 
hysteresis as shown in Fig. 121 Although relatively steep variation of the average is observed, 
the hysteresis can be seen only for the most probable values. 
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FIG. 6: The contour plot of the r-dependence of the PDFs. The thin lines represent the contour 
lines. All the contours larger than 3.43 are displayed with the same dotted lines. The bifurcation 
diagram represented with the bold solid line is added for comparison. The location of the region 
with large probability, mainly denoted with the thin solid lines, is shifted as the value of the 
parameter r is changed. However, there is no singularity on that. 
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FIG. 7: The log-log plot of the tails of the steady state PDFs for different values of r. All the tails 
are well-approximated by the same power-law. 

V. SUMMARY AND DISCUSSION 

Finally, we summarize our results and consider their implications. We applied the sta- 
tistical theory to the model of the inhomogeneous turbulent plasma where two turbulences 
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FIG. 8: Comparison of the averages with the most probable values. The averages are indicated 
with the dashed arrows and the most probable values (the peaks of the PDFs) are indicated with 
the solid arrows. The difference due to the long tail shown in Fig. [7| means that the average is not 
equal to the value which x\ takes with the largest probability. 
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FIG. 9: The r— dependence of the most probable values (solid lines). The average is also plotted 
for comparison (a dotted line). Although relatively steep variation of the average is observed, the 
hysteresis can be seen only for the most probable values. 
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well-separated in those space-time scales coexist; the "high wave number mode" turbulence 
(the CDIM micro turbulence) and the "low wave number mode" turbulence (the ITG mode 
semi-micro turbulence). We derived the development equations of the PDFs of the spectrum 
amplitudes of the electro-static potential for the characteristic wave numbers. By numeri- 
cally solving the adiabatically approximated Fokker-Planck equation, the steady state PDFs 
for the low wave number mode turbulence were obtained. Although the subcritical bifurca- 
tion is observed when the turbulent fluctuations are ignored, the shape of the PDF shows 
that the transition is smeared out by the fluctuations. It means that the approximation 
ignoring the turbulent fluctuation like the traditional transition theories could overestimate 
the range of cusp catastrophe. 

We also compared the average values with the most probable values and showed that 
these two characteristic values of stochastic nature are different due to the long power-law 
tail of the PDF. It means that the average does not mean the value which is expected to 
realize most probably and hence description of the state of the turbulent systems needs not 
only the average but also other statistical quantities like the most probable values and the 
variances. 

Consequently, these our results warn that the deterministic description of high- 
temperature plasmas cannot capture important information of the turbulent systems and 
the statistical analyses with the Langevin equations and the PDFs are inevitably needed. 
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